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We analyze dephasing in a model system where electrons tunnel sequentially through a symmetric 
interference setup consisting of two single-level quantum dots. Depending on the phase difference 
between the two tunneling paths, this may result in perfect destructive interference. However, if the 
dots are coupled to a bath, it may act as a which-way detector, leading to partial suppression of the 
phase-coherence and the reappearance of a finite tunneling current. In our approach, the tunneling 
is treated in leading order whereas coupling to the bath is kept to all orders (using P(E) theory). We 
discuss the influence of different bath spectra on the visibility of the interference pattern, including 
the distinction between "mere renormalization effects" and "true dephasing". 
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I, INTRODUCTION 

The destruction of quantum-mechanical phase coher- 
ence due to coupling of a system to an irreversible bath 
is a subject important not only because of its connec- 
tion to fundamental issues (the quantum measurement 
process and the quantum-classical transition) but also 
because of its role in the suppression of phenomena re- 
sulting from quantum interference effects, such as those 
studied in mesoscopic physics (including Aharonov-Bohm 
interference, weak localization and universal conductance 
fluctuations). Recently, the field of mesoscopic physics 
in particular has seen a revival of interest in these ques- 
tions, due to surprising experimental findings^ concern- 
ing a possible saturation of the weak-localization dephas- 
ing rate at low temperatures, that have not yet been ex- 
plained convincingly. Apart from investigations dealing 
directly with the problem of weak localization in a dis- 
ordered system of interacting electrons, several toy mod- 
els have been analyze d 2 '^ 4 'M'7ifiii? to answer the question 
whether decoherence at zero temperature is possible at 
all, contrary to the expectations based on perturbation 
theory. One of the difficulties faced by models involv- 
ing discrete levels consists in the fact that destruction of 
phase coherence for a superposition of excited states of 
finite excitation energy is perfectly possible even at zero 
temperature, due to spontaneous emission of energy into 
the bath. It is only in the zero-frequency limit of the lin- 
ear response in a system with a continuous spectrum (rel- 
evant for weak-localization and other equilibrium trans- 
port experiments) that perturbation theory suggests in 
general a vanishing dephasing rate, because then the per- 
turbation does not supply energy to the system, such that 
at T = the system is not able to leave a trace in the 
bath, which is considered to be the prerequisite for deco- 
herence. 

Some questions of interest concerning dephasing, es- 
pecially in connection with mesoscopic systems and low 
temperatures, are the following ones: How reliable is the 
simple classical picture of a phase being randomized by 
fluctuating external noised? In particular, what is the 
meaning of the zero-point fluctuations of the bath in this 




Figure 1: The double-dot "double-slit" setup, with a fixed 
phase difference ip between the two paths and under the in- 
fluence of a fluctuating environment. 



picture, as opposed to the thermal fluctuations dominat- 
ing at frequencies lower than the temperature? When 
do the former lead to "mere renormalization effects" and 
how is it possible to distinguish these from "true" de- 
phasing? Under which circumstances is the suppression 
of off-diagonal terms in the reduced system density ma- 
trix itself already a good indicator of dephasing? How 
reliable are simple arguments based on Golden Rule and 
energy conservation, related to the connection between 
dephasing and the trace left in the bath by the parti- 
cle ("which-way" detection)? When does perturbation 
theory fail qualitatively, what is the influence of non- 
Markoffian behaviour? How does the dephasing rate de- 
pend on the energy supplied by an external perturbation 
(frequencies excited in linear response, bias voltage ap- 
plied in a transport measurement)? What is the influence 
of the Pauli principle in a system of degenerate fermions? 
How strong are the qualitative differences in behaviour 
resulting from different bath spectra? 

In this work, we will present a model that is able to 
give insights into most of these questions. 

Our model represents a kind of mesoscopic double-slit 
setup. It consists of two single- level quantum dots which 
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are tunnel-coupled to two leads, with a possible phase 
difference between the two interfering paths (see Fig. . 
Due to destructive interference (at (p = ir), the tunneling 
current may be suppressed completely, provided the two 
dot-levels are degenerate and the setup is symmetric in 
the two interfering paths. Coupling the dots to a bath 
may partly destroy the phase coherence and re-enable 
the electrons to go through the setup. For a symmetric 
setup, with equal coupling strength between the bath 
and each of the two dots, mere renormalization effects 
will not be able to lift the destructive interference in this 
way. Thus, a finite tunneling current may be taken as a 
genuine sign of dephasing. This criterion for dephasing 
has been employed before in a model of dephasing due 
to spin-flip transitions in first-order tunneling transport 
through one or two dotsii, as well as for cotunneling 
through an Aharonov-Bohm ring coupled to a fluctuating 
magnetic fluj£. 

The influence of phonons on sequential tunneling 
through two quantum dots in series has been stud- 
ied experimentally in Ref. |T2- There, inelastic tran- 
sitions induced by piezoelectric coupling to acoustic 
phonons in GaAs have been essential for obtaining a 
finite current through the two off- resonant dot levels. 
This kind of setup has b een analyzed theoretically in 
Refs. llfll^lfjllfllTliallil On the other hand, we will 
be analyzing tunneling through two dots placed in paral- 
lel. Early theoretical investigations of this problem (with- 
out a fluctuating environment) include Refs. 12(11211 Re- 
cently, a parallel-dot tunneling setup has been realized 
experimentally in Ref. HE with an emphasis on spec- 
troscopy of the "molecular states" of the double-dot sys- 
tem (with inter-dot tunneling present). In our model 
of an interference setup, we choose to describe a sit- 
uation without tunneling between the dots (but with 
Coulomb-repulsion). In addition, we want to concen- 
trate on interference effects in the orbital motion and 
therefore consider the case of spin-polarized transport. 
This model - in the absence of a fluctuating environment 
- has been investigated previously in Ref. Other 
recent theoretical works concerning tunneling through 
dots in a parallel geometry have mostly investigated spin 
and Kondo physic g 24 ' 2 §i 2 fi , but also dephasing by spin- 
flip transitions^. Some works have treated the influ- 
ence of phonons in tunneling interference structures 2 ^ 2 ^, 
but no systematic discussion of dephasing and the visi- 
bility of the interference pattern has been given. Some 
while ago, dephasing by nonequilibrium current noise has 
been investigated experimentally^ and theoretically^ in 
a setup with a single quantum-dot placed into one arm 
of an Aharonov-Bohm interferometer. 

Our analysis of dephasing in sequential tunneling 
through a double-dot will take into account the system- 
bath coupling exactly, while we treat the tunnel-coupling 
only in leading order. The presence of the Fermi sea in 
the leads introduces some aspects related to the Pauli 
principle and to the behaviour of systems with a contin- 
uous spectrum that cannot be analyzed in simpler models 



of dephasing in discrete systems coupled to a bath. 

The work is organized as follows: After setting up the 
model (Sec. [H}, we will present a qualitative discussion 
of its main features ijllljl . In particular, we will discuss 
the relation between entanglement, dephasing and renor- 
malization effects. Subsequently, we derive a general for- 
mula for the tunneling decay rate of an electron that has 
been placed on the two dots in a symmetric superposi- 
tion of states (Sec. IIVII . This is done by building on 
the concepts of the P(E) theory of tunneling in a dissi- 
pative environmen t^ 1 "^ 2 . Following this, we will evaluate 
the dependence of the tunneling rate on the bias voltage 
and the bath spectra (Sec. EJ. We will interpret the 
results in terms of "renormalization effects" and "true de- 
phasing" (Sec. IVIji . Building on these sections, we will 
finally derive a master equation for the case of weak tun- 
nel coupling (Sec. IVIIj) . which allows us to calculate the 
sequential tunneling current as a function of bias voltage, 
temperature, and phase difference (Sec. IVIIIl! . 

The most important results derived in this work are 
the following: Equation 111 -ill is the general expression for 
the phase-dependent tunneling decay rate in presence of 
the fluctuating environment. It forms the basic input 
for the master equation (Eqs. l(56jl - (l58ll % ). that describes 
sequential tunneling through the double-dot, where the 
resulting current can be obtained from Eq. IjfiOJI . The vis- 
ibility of the interference pattern, which is defined by the 
phase-dependence of the current, is given in Eq. It 
is connected with the visibility obtained from the phase- 
dependence of the tunneling rate itself (Eqs. ill fit . j!7t ). 

II. THE MODEL 

We consider a Hamiltonian describing two degenerate 
single-level quantum dots, with respective single-particle 
states |+) and |— ) (spin is excluded for simplicity, since 
we are interested in dephasing of the electronic motion) . 
Each of them is tunnel-coupled to two electrodes (with 
the same strength for both dots) , but involving a possible 
phase difference between the tunnel amplitudes (see Fig. 
nj. In addition, the potential difference between the two 
dots is given by a fluctuating field F, whose dynamics is 
derived from a linear bath. It represents the fluctuations 
due to phonons or Nyquist noise. The system-bath cou- 
pling strength is taken to be the same for both dots, while 
the sign is opposite, such that the bath can distinguish 
between an electron being on dot |+) or |— ): 

H = e{h + + n_) + F(h + — ri_) + Un + n_ + 

H L + Hr + Hb + V (1) 

Here n± are the particle numbers on the two dots 
(equal to or 1). The bath Hamiltonian Hb describes a 
set of uncoupled harmonic oscillators. It governs the dy- 
namics of the fluctuating potential F, which is assumed 
to be linear in the oscillator coordinates. The coupling 
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between electron and bath is of the form of the "indepen- 
dent boson modeF'^i. For the case of exactly one electron 
on the double-dot, and in the absence of tunneling, it cor- 
responds to a spin-boson model with "diagonal coupling". 
In this model, no transition between different levels is 
brought about by the bath, such that pure dephasing re- 
sults. U denotes the Coulomb repulsion energy, which 
we will take to be so large that double-occupancy is for- 
bidden. Note that the degeneracy of the two dot-levels is 
important in the following: It is necessary to ensure com- 
plete destructive interference at <p = 7T (compare also the 
discussion in Sec. IVIIj) . 

The terms H j, and Hr contain the energies of the elec- 
trons in the left and right reservoirs: 

Hl(r) = £k& L(R)k d L(R)k ■ (2) 
k 

The tunneling between the dots and the leads is described 
by V = Vl + Vr, with 

V R = Y, + + h.c. (3) 

k 

for the right junction, and 

VL=Y t k & Lk(d+ + d-) + h.c. (4) 

k 

for the left junction. 

Here d± are the annihilation operators for the two dots 
(n± = d\d±) and the phase-factor of e l(p controls the in- 
terference between tunneling events along either the up- 
per or lower path. The tunneling phase difference might 
be thought of as arising due to the Aharonov-Bohm phase 
from a magnetic flux penetrating the region between the 
quantum dots. 

Note that the tunneling matrix elements t k are as- 
sumed not to depend on the dot state |+) or |— ) in our 
model. This means that the dots are close enough such 
that they couple to the same point on the lead electrodes, 
to within less than a Fermi wavelength. Obviously there 
could be no appreciable interference effect if the dots were 
separated by some larger distance (in which case the k- 
dependence of matrix elements would be different for the 
two states). The same idealized ass umption underlies 
several similar models (see, e.g., Refs. Illl23l2fil) . The ef- 
fect of an arbitrary dot separation has been discussed in 
some detail in Ref. [2lJ. 

The present model, without the bath, has been ana- 
lyzed previously in Ref.E^l (see also Sec. IV. C of Ref.0)- 
There, an orbital type of Kondo effect was found in equi- 
librium, for ip = 7r, when the level energy was below the 
chemical potential. This arises because at tp = ir there 
are two states of the double-dot that couple only to the 
left and the right lead, respectively (denoted by |e) and 
\o) in the following). These degenerate states form the 



pseudospin responsible for the Kondo effect. However, 
that mechanism will be irrelevant for our analysis, as 
we consider the transport situation where the (renormal- 
ized) level energy lies between the chemical potentials of 
the left and the right lead. Therefore, the degeneracy is 
effectively lifted by the bias voltage (which will be as- 
sumed to be much larger than the tunneling rate), and 
only the state coupling to the left lead would be occupied 
at tp = 7T. 

III. QUALITATIVE DISCUSSION 

In this and the following three sections, we first ana- 
lyze the escape of a single electron into the right lead, 
where the electron is assumed to start out in a symmet- 
ric superposition of the two dot levels, which has been 
formed by an electron tunneling onto the dots from the 
left lead. In the situation without any bath, this is the 
state |e) = (|+) + \-))/V2. 

Without dephasing, the tunneling decay out of state 
|e) is made impossible in the case of perfect destructive 
interference at tp = 7T, while maximal constructive inter- 
ference is present for tp = 0. It should be noted that the 
attribution of the phase factor to one of the tunnel cou- 
plings represents a certain choice of gauge, which affects 
the wave functions in the following discussion but none 
of the physically observable quantities that are derived 
as a result of the master equation in Section IVIII 

For simplicity, we will assume a zero-temperature sit- 
uation throughout the following qualitative discussion, 
with a bias eV > applied between the two dots and the 
lead in such a way that the electron is allowed to tun- 
nel into the lead (see Fig. [21 . In addition, since we will 
describe the tunneling decay to the right, we will only 
consider the coupling Vr to the right lead in this section 
and drop the index R for now. 

Without the bath and for perfect constructive interfer- 
ence {tp = 0), the tunneling decay rate T will take on its 
maximum value of 2Fo , with 

r = 2ttZ^ (|t fc | 2 ) , (5) 

where D is the lead density of states at the Fermi en- 
ergy, (\tk\ 2 ^ is the angular average of |tfc| 2 at this energy. 

The bias voltage V does not enter in this case, as long as 
it is positive (permitting decay). For tp = n, T vanishes 
due to perfect destructive interference. In general, we 
have: 

r = r (i + cos^). (6) 

If the bath is included in the description, the following 
happens: 

First of all, the energy of a single extra electron on any 
of the two dots will be renormalized from its initial value 



4 




Figure 2: The ground state |x+) (Ix-)) which the bath as- 
sumes in the presence of an electron on dot |+) ([—)), shown 
schematically for a single bath oscillator (see main text). Af- 
ter the electron has tunneled into the lead, |x~) becomes a 
superposition of excited states (dashed), while the state |xo) 
represents the ground state of the bath in the new potential. 



of e, since the bath relaxes to a ground state of lower en- 
ergy in presence of the electron. We will assume that the 
value of e has been chosen exactly to compensate for this 

energy change, which is given by — J °° du; (Pf\ l u 

(see App. EJ. Then, the energy of an electron on the 
dot (and the bath in its new ground state) is the same 
as that of the electron being in the lead, at the Fermi 
energy of = (for V = 0). 

Tunneling of an electron from the dots to the lead will 
not change the bath state, but it will displace the origin 
of the harmonic oscillators comprising the bath, since 
the coupling to F is switched off (n+ — h- changes to 
zero). Therefore, the original ground state of the bath 
(in presence of the electron) will become a superposition 
of excited states in the new bath potential (in absence 
of the electron; see Fig. [2j. On the other hand, since 
energy conservation has to be fulfilled with respect to 
the total energy of the electrons and the bath before 
and after the tunneling event, only those excited bath 
states can be reached whose energies are not greater than 
eV, the energy supplied to the electron by the bias volt- 
age. This leads to the Coulomb-blockade type suppres- 
sion of the tunneling rate at low bias voltages, for <p = 0. 
Physically, this effect is just the same as that described 
by Franck-Condon overlap integrals evaluated between 
vibronic states for electronic transitions in molecules. 
Qualitatively, this effect is independent of the interfer- 
ence setup, since it already occurs for tunneling through 
a single dot coupled to a bath. 

In contrast, for the case of destructive interference 
(ip = 7r), the bath may actually enhance the tunneling 
rate from its initial value of 0, since it partly destroys 
the phase coherence that is a presupposition for perfect 
interference. An electron coming from the left lead will 
form the following entangled state with the bath, instead 
of the symmetric superposition |e) = (|+) + \ —))/y/2: 



(|+>|x+> + |->IX-»/>/2. (7) 

Here the states |x±) denote the respective ground 
states of the bath for a bath Hamiltonian given by 
Hb ± F, which are related to each other by a parity 
transformation (This also means we assume by defini- 
tion there to be no phase factor between these states; 
e.g. both may be assumed to have real-valued positive 
wave functions) . Actually, the entangled state considered 
here will be formed only if the electron is given barely 
enough energy to enter the double-dot at all (i.e. chem- 
ical potential of the left lead infinitesimally larger than 
the renormalized level position). Otherwise, excited bath 
states may be created even at this step. These complica- 
tions will be taken care of in the complete discussion of 
the sequential tunneling current (Section IVIIJl . There, it 
will turn out that the tunneling decay rate derived in the 
following, based on our physically motivated ansatz Q, 
is exactly the rate that enters the full master equation. 
Thus, we proceed with the ansatz Q for the initial en- 
tangled state, in order to calculate the rate for such an 
electron to tunnel into the right lead. 

The bath measures (to some extent) which dot the 
electron resides on, such that the reduced system den- 
sity matrix (for the electron on the two dots) becomes 
mixed and its off-diagonal elements get suppressed by 
the overlap factor (x+lx-)- Put differently, the phase 
factor between the two dot states in the wave function 
of the electron (initially equal to +1) becomes uncertain. 
Therefore, there is a finite probability of 

Po = (l-<X+|X-»/2 (8) 

to find the electron in the antisymmetric (odd) state 
\o) = (|+) — |— ))/v / 2- At <p — 7r, where tunneling de- 
cay of the symmetric superposition |e) is blocked due to 
destructive interference, the state \o) is allowed to decay 
into the lead, at the maximal rate of 2To. In this way, 
the interference-induced blockade of electron tunneling is 
lifted by dephasing. 

However, this simple picture is true only for large 
bias voltages, when energy conservation permits any fi- 
nal state of the bath after the tunneling event. If the 
maximum energy supplied to the electron is limited, the 
suppression discussed above (for the case of ip — 0) will 
apply again. In particular, if the bias voltage is turned 
to zero, energy conservation only allows the state |xo) to 
be reached, which is the ground state of the bath in the 
absence of any electrons on the dots. Then, the tunnel- 
ing rate is exactly zero again, despite the fact that the 
reduced density matrix of the electron may be mixed to 
a strong extent. The reason is the following: When the 
overlap of the entangled state with the state |xo) is 
taken, the two overlap factors (xolx+) an d (xolx-) turn 
out to be the same, if the coupling of the bath to the 
two dots is symmetric (i.e. of equal strength, only of 
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opposite sign), which we have assumed in writing down 
the Hamiltonian, Eq. Q. Therefore, the electronic state 
resulting from the projection of @ onto |xo) is equal 
to the symmetric combination, whose decay is forbidden. 
Thus, the combination of energy conservation and Pauli 
blocking prevents a finite tunneling rate at zero bias volt- 
age, in spite of the mixed state of the electron coupled 
to the bath. In this limit the entanglement between elec- 
tron and bath only leads to renormalization effects (such 
as the change in tunneling rate), but not to genuine de- 
phasing. If the coupling were asymmetric, then destruc- 
tive interference could be lost even without dephasing 
(merely due to renormalization), just as it would be the 
case for initially asymmetric bare tunnel couplings. That 
is why the asymmetric case is uninteresting for our pur- 
poses of distinguishing renormalization effects from real 
dephasing. 

However, whether we are indeed able to claim that de- 
phasing actually vanishes in the limit of low bias voltages 
will depend on the behaviour of the tunneling rate as a 
function of V and on the comparison of the cases <p = 
and ip = ir. Here, the bath spectrum, and, above all, 
its low-frequency properties, enter. In order to be able 
to discuss T(V, ip) quantitatively, we will make use of the 
concepts of the P(E) theory of tunneling in a dissipative 
environment. 



IV. DECAY RATE AND CONNECTION TO 
P{E) THEORY 

The tunneling rate T will be calculated using the stan- 
dard Fermi Golden Rule, i.e. lowest order perturbation 
theory in the bare tunneling rate To, but taking into ac- 
count exactly the bath coupling. In deriving the formula 
for r, it turns out to be useful to assume that the bath 
oscillators do not get shifted in the tunneling event (un- 
like the qualitative considerations from above), but it is 
rather the bath states which get displaced (in the op- 
posite direction). Obviously, this amounts to the same, 
as long as we are interested only in overlap integrals of 
different bath states after the event. To that end, we in- 
troduce the displacement operator exp(i<f)), which trans- 
forms the bath ground state of Hb into that of Hb + F. 
Here <j> is a suitable hermitian operator that is linear in 
the bosonic variables of the bath. In fact, this amounts to 
performing the canonical transformation of the indepen- 
dent boson model22, see AppendixEl In terms of the two 
dot states + and — , we have F + — F and F_ = — F, as 
well as 4>+ — <f> and (j>- = — </>. The transformation elim- 
inates the system-bath coupling from the Hamiltonian, 
but gives rise to modified dot operators d' ± = e ±l ^d± in 
the transformed tunnel Hamiltonian V' R (see Eq. QA7j) l. 

We will assume the tunnel-coupling to be sufficiently 
weak, such that we can use lowest-order perturbation the- 
ory to calculate the tunneling decay rate: 



r = 27r]T|(/|V^) S(E f -Ei), 



f 



(9) 



where the initial state \i) is given by the configuration 
involving the electron residing in the symmetric superpo- 
sition on the dots, the unperturbed Fermi sea in the lead 
and the bath in its ground state \ib)- The bath ground 
state has become independent of the position of the elec- 
tron, due to the above-mentioned transformation. At 
finite temperatures, an additional thermal average over 
the initial bath state and the initial state of the electrons 
in the lead has to be performed. The energies and eigen- 
states refer to the Hamiltonian without tunnel coupling. 
Applying the new tunneling Hamiltonian V' R to the initial 
state, we obtain the following expression: 

r = 7r^ \t k \ 2 (l- f(e k + eV))x 

kJ B 

f B \e+^ + e^e-^\ lB )\ 2 5{Ef - E? + e k ) , (10) 



Here /(•) is the Fermi function (for chemical potential 
equal to zero), and Ef fi are the energies of the initial 
and final bath states. The energy supplied to the bath 
is equal to the energy lost by the electron (given by — 
since the renormalized dot energy is zero) . Following the 
usual derivation of the P{E) theor yi? 1 ^ 2 , we express the 
energy-conserving 5 function as an integral over time and 
also replace the sum over lead states k by an integral over 
the energy E = —e k supplied to the bath, finally yielding: 



-e iEt x 



r = r / dE(l-f(-E + eV)) ^< 

J — oo J —oo ^ 

'( e -#(t) + e -^e^ (4) )(e^ + e iip e~ trp ) ) (11) 



' — OO 

1 



For the case of arbitrary temperature, the brackets de- 
note a thermal average over the initial bath state \ib)- 
We introduce the definitions: 



JEt „±(d>(t) 



(12) 



This permits us to write down our final result for the 
tunneling decay rate in terms of P(_)(E): 



r = r / dE(l-f(-E+eV))(P{E)+cos(ip)P-(E)) 

J —oo 

(13) 

The formula given here constitutes the basic expres- 
sion for the decay rate as a function of bias voltage and 
interference phase ip. It represents the appropriate mod- 
ification of Eq. 10 in presence of a bath. 
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Note that for the slightly more general case of arbitrar- 
ily correlated fluctuating potentials F + and F_ attached 
to the dots (i.e. an interaction of the form F + n + +F_n_), 
the function P- (E) would contain the cross-correlator of 
the associated phases <f>+ and 0_ , while P(E) would de- 
pend on the autocorrelator of </>+ or <fi- (assumed to be 
the same, for the setup to remain symmetric). In con- 
trast to the model treated here, such an interaction would 
also involve fluctuations of the sum of energies of the dot- 
levels. However, they would only add to the renormaliza- 
tion effects mentioned previously and do not contribute 
to dephasing by themselves, since such fluctuations can- 
not distinguish between the two interfering paths. 

By using the definitions 

7(-) = T JdE{\- f(-E + eV)) P ( _j (E) , (14) 
we can write 



r = 7 + cos(y)7_ 



(15) 



The strength of the dependence of T on the phase <p 
may be taken as a signature of phase coherence in our 
model. We define the "visibility" of the interference pat- 
tern in the usual way, by 



^ — (J^jnax ^min ) / (J^max ~t~ ^min)i (-^) 

which is equal to the ratio 



■1_ 

1 



(17) 



The visibility v will be 1 whenever the destructive in- 
terference is perfect, and it is zero if there is no depen- 
dence of r on ip. 

The effects of the bath on the decay rate are encoded in 
the functions P{E) and P-(E), whose general properties 
we will discuss now. In the next section, we will evaluate 
them for different types of bath spectra. 

As usual, the function P{E) describes the probability 
(density) that an electron will emit the energy E into the 
bath while tunneling into the lead. It is real, nonnegative 
and normalized to unitj«2i*22. 

At large times \t\ — > oo, the correlation function 

(3>(t)4>) in the exponent of the integral l(T2|l will de- 
cay towards zero, for a continuous bath spectrum. This 
means that the integrand of P(E) approaches the value 

of z = cxp(— ^0 2 ^}), starting from 1 at t = 0, There- 
fore, P(E) contains a "quasiparticle 5 peak" of strength 
z at E = 0, if z does not vanish. It corresponds to the 
probability z of having no energy transfer at all from the 
electron to the bath (similar to the recoil-free emission of 
a 7 ray by a nucleus inside a crystal, i.e. the Mossbauer 
effect). 



The function P-(E) in front of the cos(^) term in Eq. 
l(T3l) is different: The integrand of P-(E) will increase 
at large times, towards the value of z, starting from z 2 
at t = 0, The function P-{E) is real- valued (because 

of (^cj)(t)<fij = /(^>(i)\ ), but it can become negative. 

Therefore, it cannot be interpreted as a probability den- 
sity, in contrast to P(E). Its normalization is given by: 



J dEP_(E) = 



(18) 



If z is nonzero, P-(E) also has a S peak at E = Q, of 
weight z, just as P(E). As a consequence, in the case of 
destructive interference (tp = %), the tunneling rate T at 
V —>■ 0, T = still vanishes even in the presence of the 
bath, since the S peaks contained in P(E) and P-(E) 
cancel exactly in the integral l(T3)l . The physical reason 
for this coherence has been discussed at the end of the 
previous section. 

In the case of constructive interference (ip = 0), at T = 
and for V — ► 0, the integration over E will only capture 
the 5 peaks contained in P^(E), yielding T = 2zT n . 
Thus, the tunneling rate is suppressed by the constant 
factor z from its noninteracting value. However, this may 
be interpreted as a mere renormalization of the effective 
tunnel coupling, since the visibility v of the interference 
pattern is still equal to unity. In order to connect this 
result to the qualitative discussion from above, we note 
that the overlap of the two different bath ground states 
that are adapted to the absence or presence of an electron 
on dot ±, is given by: 



(Xo|x±) = (xo 



Xo 



exp(-(0 2 )/2) = z 



,1/2 



(19) 

Therefore, the magnitude squared of this overlap, that 
determines the probability of tunneling without exciting 
any bath mode, is equal to z. 

On the other hand, for sufficiently large bias voltages 
(much larger than the cutoff frequency of the bath spec- 
trum), the normalization conditions for Pr^[E) yield 



r (l + z 2 cos(<^)) . 



(20) 



The visibility is given by v 



In this limiting case, 



where the restrictions due to energy conservation and the 
Pauli principle are no longer important, the tunneling 
rate T at the point ip = ir of destructive interference 
does not vanish. It takes the value To(l — z 2 ), which is 
small if the effects of the bath are weak (z near to 1) and 
is equal to one half the ideal maximum value 2]?o for a 
bath that is sufficiently strong to destroy phase coherence 
completely (z = 0), leading to an incoherent mixture of 
symmetric and antisymmetric states on the two dots. In 
the latter case, the visibility vanishes (even for arbitrary 
voltages), since then P-{E) is equal to zero, which makes 
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r independent of p. This will be true for the Ohmic bath, 
to be discussed in the next section. 

As explained above, the reduced density matrix of the 
electron on the dots coupled to the bath predicts a finite 
probability of P a = (1 — (x + |x_))/2 to find the electron 
in the antisymmetric state if one starts out from the sym- 
metric superposition before coupling it to the bath. The 
overlap factor of the bath states involved in this proba- 
bility can be expressed as 

<X+IX-) = (xo|(e-^) 2 |xo)^ 2 . (21) 

Comparing with the result T((f = tt) = r (l — z 2 ) given 
above, it may be observed that the decay rate at suffi- 
ciently large bias voltages is indeed determined directly 
by the probability to find the electron in the state whose 
decay is not forbidden by destructive interference (as has 
been argued already at the end of the previous section, 
near Eq. (jHJ). It is only in this limiting case, where 
an arbitrary amount of energy is available for excitation 
of the bath, that the suppression of interference effects 
in the transport situation is correctly deduced from the 
electron's reduced density matrix in the presence of the 
bath. Formally, this holds because the sum over final 
bath states fs in Eq. itTTHl is not restricted any more and 
corresponds to the insertion of a complete set of basis 
states. Thus, one obtains, directly from Eq. 

T = ^-( X+ + e-^x-\x+ + e iv X-) , (22) 

which reduces to Eq. 12011 when the overlaps are eval- 
uated, using Eq. 12 111 . Physically, the case of high bias 
voltage corresponds to a kind of infinitely fast von Neu- 
mann projection measurement that determines the state 
of the electron, revealing the fluctuations due to the bath. 
In contrast, at low bias voltages (low energy supply), a 
kind of "weak" measurement is carried out that takes a 
longer amount of time, such that only the low-frequency 
fluctuations of the bath are important for dephasing. 



V. EVALUATION FOR DIFFERENT BATH 
SPECTRA 

We will restrict the discussion to T = at first. 

The simplest example for the bath is a single harmonic 
oscillator of frequency oj. This offers an approximate de- 
scription of the interaction with optical phonon modes 
("Einstein model"). In this case, P(E) and P-{E) can be 
obtained easily by expanding the exponential in a Taylor 

series and using (^$(t)^J — (4> 2 ^ ex P(~ iut), before the 

integration over time is performed. For P(E), the result- 
ing series of 5 peaks at harmonics of oj corresponds to 
all possible processes where the electron emits any num- 
ber n of phonons into the bath while tunneling into the 



lead. The expression for P- (E) is the same, apart from 
alternating signs in front of the 5 functions: 

P { _ ) (E) = zJ2^r-S(E-nu;). (23) 

Thus, every process involving the transfer of an even 
number of quanta to the bath will not ruin the destructive 
interference at <p = 7T, since the corresponding contribu- 
tions from P{E) and P-{E) cancel in Eq. lfTH)i. This 
is because the coupling between electron and bath is of 
the type (n+ — h^)F , which gives a different sign of the 
interaction amplitude for a phonon emission process, de- 
pending on the dot. Therefore, the amplitude of emission 
of an even number of phonons will not depend on the dot, 
it is insensitive to the state of the electron, and the am- 
plitudes of the electron tunneling from |+) and |— ) will 
still interfere destructively. 

In contrast, emission processes involving an odd num- 
ber of quanta introduce a negative sign for an electron 
starting in state |— ), "detecting" the path (or rather, the 
initial state) of the electron and interfering constructively 
with the processes from |+). This lifts the destructive in- 
terference and makes r ^ at <p = 7r. However, below 
the frequency oj of the oscillator, destructive interference 
at tp = ir is still perfect since no quantum can be emitted, 
while the magnitude of T at p — is renormalized by the 
factor z, as has been discussed above in general for the 
limiting case V — ► 0. The same holds true for any bath 
with a finite excitation gap, at T = 0. This is shown in 
Figs. 31 and to be discussed in the next section. 

We now pass on to arbitrary bath spectra. At first, 
we will cover the case z / ("weak baths"), when we 
can apply perturbation theory to discuss the behaviour 
of P(_)(i5) at low energy transfers E (and, consequently, 
that of T at low voltages). A Taylor-expansion of the 
exponent in Eq. ((12)) yields: 

Z7r „=o n - 

= ^E^f((#) *•••*(#) )(^) ( 24 ) 

The repeated convolution product contains n times the 

correlator ($4^ , for n = it is to equal S(E), and the 

negative sign holds for P_ (E) . 

For the following discussion, we prescribe the spectrum 
of the fluctuating potential F to be a power-law in fre- 
quency oj (at T = 0), with exponent s: 

(PF^' = 2aoj c (jpj 0{oj c -oj)6{oj), (25) 

The dimensionless parameter a characterizes the bath 
strength. In order to be able to rely on perturbation 
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theory, we have to ensure z > 0. Since 

^PP\ /uj 2 , the variance of the fluctuating phase, 

will be finite only for s > 1 (at T = 0, otherwise s > 2). 
In that case, we have z — exp(— 2a/ (s— 1)). This means 
the perturbative analysis presented above is restricted to 
a super-Ohmic bath, s > 1. The case of the Ohmic bath 
will be discussed separately further below. 

After keeping only terms up to second order in the 
expansion of Pr_\(E) given in Eq. 124L we get 



P{E)+P^E)=z{25{E) + ({^4>) *(#) ){E) + ...), 

(26) 

for the symmetric combination, that will determine the 
prefactor of 1 + cos(</j) in the expression for T, Eq. ill .'"it . 
and 



P{E) - P-{E) = 2z 



(27) 



for the antisymmetric combination (determining the 
prefactor of 1 — cos(<^)). Inserting these into 11 .'111 , using 

the power law for (^><p) = (pF^j /uj 2 given by p5f . 

and performing the energy integrals, we find, for suffi- 
ciently low voltages (2a(eV/w c ) s <C 



To r , , a 2 C s feV 

r » ^z{ i + cos ^ i + - — s - — 

2 {s -1) \u c 



1): 



2(6-1) 



(1 - cos(y>)) 



2a /eV 



s — 1 \ UJ, 



}.(28) 



The numerical prefactor C' s is defined as Jy(y(l — 
y)Y- 2 dy. 

From Eq. H28I1 , we see that the destructive interference 
at cp = 7T is perfect at V = 0, but gets lifted when increas- 
ing the bias voltage, with a power V s ' 1 . In contrast, the 
decay rate F at <p = starts out from the constant value 
of 2zTq and grows as V r2 ( s ~ 1 - ) . Therefore, the visibility v 
starts out at 1 for V — but decreases as: 



1 - 



4a / eV 



s-1 



(29) 



For s I 1, the range in bias voltage V where these ap- 
proximate expressions hold shrinks to zero (at constant 
a and uj c ). At s = 1, i.e. for the Ohmic bath, the prob- 
ability z of not emitting energy into the bath vanishes 
completely. As discussed above, this means that there 
is no (/^-dependence at all in T, and, consequently, the 
visibility is zero at all bias voltages. Furthermore, the 
tunneling rate vanishes for eV — ► 0, even at tp = 0. 
This is the well-known Coulomb-blockade type of be- 
haviour for tunneling in presence of Ohmic dissipation^!. 
At higher bias voltages, the blockade is removed and T 



grows towards Tq. The growth at low voltages is deter- 
mined by the power-law behaviour of P(E), which rises 
where the exponent is determined by the 



2a rp2a-l 



as CUJ C " J £ 

bath-strength rather than the exponent s — 1 of the bath 
spectrum. The dimensionless prefactor c must be found 
from the normalization condition for P(E) and depends 
only on a (and the type of cutoff in the bath spectrum) . 
Therefore, in the case of the Ohmic bath we have, at low 
V and T = 0: 



2a \ uj r 



2a 



(30) 



Finally, we briefly discuss the case of finite tempera- 
tures, T > 0. 

In that case, the variance of (f> is given by 



which yields 

(T=0) 



(T=0) 



coth 



(-) 

\2TJ 



(31) 



eV - 1 



dy. (32) 



The approximation of extending the integral to infinity 
holds for temperatures much smaller than the bath cut- 
off uj c . This formula gives the temperature-dependence 

of the renormalization factor z = exp (^- ($ 2 j\. The 

second integral diverges for s < 2, because z = for 
these cases, in contrast to T — where 2 = only for 
s < 1. Again, this results in complete absence of the 
interference effect in the tunneling rate T(V, <p) (because 
P-{E) vanishes). It may seem surprising that an in- 
finitesimally small temperature can yield such a drastic 
qualitative change (for 1 < s < 2), compared to the 
zero-temperature case, since the difference should be ob- 
servable only at very large times t jp l/T. However, it 
must be remembered that our analysis is carried out for 
the limit To — * 0, where the average decay time of the 
given state is inifinitely large. In other words, the limits 
T — > and ro — > do not commute for such relatively 
strong baths. At finite To, the transition from one to the 
other case should turn out to be smooth, but this goes 
beyond the present analysis. 

Apart from the change in z with temperature, there 
are two other important differences to the case T = 0: 
First of all, even at V — > the electron may emit energy 
into the bath, due to the thermal smearing of the Fermi 
surface in the lead (lifting of Pauli blocking). Secondly, 
it may now also absorb some energy during the tunneling 
process. Both facts will, in general, lead to a finite tun- 
neling decay rate at tp — n, V — ► for any bath, where, 
at T — 0, the rate had vanished in any case. 

We can approximate the visibility v at V — > and fi- 
nite T by using the expansion ll24^l . Inserting the result- 
ing expressions for 7(_) 111 111 into v = 7-/7, we obtain 
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v(T,V->0)nl-4[de(M)j(e). (33) 

We evaluate the integral for a power-law bath spectrum 
in the limit T -C uj c : 



de (M) m 



2aui 



This yields: 



de 



T=0 



o sinh(/?e) 

fS -2 

de- 



o sinh(/3e) 



(34) 



A_ 




l-v(T,V->0) w 32a — 



M (i-2-)r(*-i)c(*-i) ; 



(35) 

where T is the Euler gamma function, and £ the Rie- 
mann zeta function. Therefore, the decrease of the vis- 
ibility with increasing temperature T (and V — * 0) is 
governed by the same power-law as that for increasing 
bias voltage V at T = 0, see Eq. (j29]) . 



Figure 3: The bath spectrum \PP/ (bottom) and the re- 
sulting functions P(E) (top) and P-(E) (middle), plotted vs. 
energy E, for different baths. Energies are measured in units 
of the "bath cutoff" u) c . Energy axis is the same in all panels 
(starting at E = 0, horizontal tick distance: 1); vertical tick 
distance in all panels is 0.5. a: s = 1.5, a = 0.25; b: "acoustic 
phonons", s = 3, a = 1; c: "optical phonons", Bath with gap; 
d: s = 1, a = 0.25; e: s = 1, a = 0.75 (d,e are "Ohmic" baths 
of different strength, z = 0) 



VI. DISCUSSION OF THE RESULTS 

The following discussion relates to the results obtained 
for T = 0, that are plotted in the figures. 

In Fig. several different types of bath spectra 

FFj are shown. Cases (a),(b),(d) and (e) are power- 

laws of the form given in Eq. ll25)l . for a cutoff fre- 
quency of uj c — 1. The last two (d,e) are of Ohmic 
type (s = 1, 2 = 0), which corresponds physically to 
gate voltage fluctuations due to Nyquist noise. Case 
(c) represents a bath with an excitation gap (for exam- 
ple optical phonons), with a spectrum given by an in- 
verted parabola. In the limit of infinitely small spectral 
bandwidth, it would correspond to the single harmonic 
oscillator (Einstein mode) discussed above. Case (b), 
with a bath spectrum rising as w 3 , corresponds to the 
experimentally relevant case of piezoelectric coupling to 
acoustic phonons, which was determined to be the ma- 
jor inelastic mechanism in the experiments of Ref. ^3 on 
double-dots in GaAs (see Ref. [la for a theoretical anal- 
ysis deriving this spectrum for wavelengths larger than 
the dot distance). The spectra for the first three cases 
(a,b,c) have been chosen to give the same renormaliza- 
tion factor, z = 1/e. The same figure shows the resulting 
functions P(E) and P-{E). These have been obtained 
using the integral equation described in Refs. I3HI35I We 
recall that the low-energy behaviour of P(E) is given by 

= (FF) E / E2 for the cases with z^0, where 



// 
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Figure 4: Decay rate F as a function of bias voltage V for the 
case of constructive interference (^3 = 0), at T = 0. Curves 
correspond to different bath spectra shown in Fig. Dashed 
lines correspond to approximation Eq. 1281 . The initial 
Coulomb-blockade type suppression to a value of F/2ro = z 
(z = for the Ohmic bath d,e) is lifted with increasing bias 
voltage, saturating at r/2r = (1 + z 2 )/2. Inset depicts en- 
ergy diagram with definition of bias voltage for this situation. 
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Figure 5: Decay rate F(V) for the case of destructive interfer- 
ence ((p = 7r), at T = 0. Dashed lines refer to Eq. 1281 . Due 
to dephasing, the decay rate becomes finite at finite voltages, 
saturating at r/2r = (1 - z 2 )/2. For the Ohmic bath (d,e) 
the dependence is exactly equal to that for ip = (Fig. ^J. 




eV/(o c 



Figure 6: Visibility v = (F v= o — Y v=v )/{T v= o + T v=7r ) as a 
function of bias voltage V for different bath spectra (see Fig. 
EJl. For the Ohmic bath (cases d,e) v = 0. Dashed lines cor- 
respond to Eq. 1291 . Inset illustrates change in interference 
pattern F((p) upon switching on the interaction with the bath. 



perturbation theory may be applied. In case (c), the al- 
ternating signs of the different contributions to P-(E) 
may be observed, whose physical meaning has been ex- 
plained above for the limiting case of the harmonic oscil- 
lator. 

We now briefly mention some numerical estimates for 
the bath strengths as they may occur in experimental 



situations. 

In GaAs, the lack of inversion symmetry leads to 
piezoelectric fields proportional to the lattice deforma- 
tion, whose effect on electrons at low frequencies is much 
stronger than that of the usual deformation potential 
(where it is only the potential that is proportional to 
the deformation). For the piezoelectric coupling^ to 
acoustic phonons in GaAs, one finds (compare Ref. Hfl) 
/ - -\ T=0 

(^(/)(f)j = W to / (c s / d) 2 for uj <C Cs/d, where c s « 

5 • 10 3 m/s is an estimate for the average velocity of 
longitudinal sound waves in GaAs, and d denotes the 
distance between the quantum dots. We obtain W — 
const- (ehu/iTr) 2 /(hpcl), where ehu = lAeV/nm is the 
single piezo-electric modulus in the cubic Td structure of 
GaAs and p = 5.3 • I0 3 kg/m 3 the mass density. The 
numerical constant is of order 1 and accounts for the de- 
tails of the sound wave dispersion relation as well as the 
orientation of the crystal axes with respect to the vector 
separating the quantum dots. Inserting these values, W 
is found to be on the order of 0.01. In order to obtain 
the renormalization factor z, the spectrum (4>4>^j must 

be integrated over all frequencies (see above), i.e. up 
to the cutoff frequency w c - The effective cutoff frequency 
uj c cc c s /do is determined by the extent do of the dot wave 
functions (for e?o = lOOnrn one obtains w c ~ 50 GHz). 
Given the present values, and assuming c?o ~ d, this leads 
to estimates for J (^4>^j duo on the order of 0.01, yield- 
ing z = cxp(— ^ 2 ^) near 1. Note that the distance d 
between the dots cancels in the estimate for z, as long as 
the cutoff frequency is assumed to be given by uj c oc c s /d. 
However, as uj c might be considerably larger than c s /d 
(if do <C d) , one could also obtain a z that deviates more 
strongly from unity. 

For the Ohmic bath, we may imagine the quantum 
dots placed inside a capacitor C connected to a circuit 
of resistance R, such that the potential difference 2F be- 
tween the dots would be given by the fluctuating volt- 
age drop across the capacitor. This leads to a bath 

spectrum ^FF^ = n(R/R Q )h 2 uj/{l + {RCuo) 2 ), with 

Rq = h/e 2 the quantum of resistance. Therefore, the di- 
mensionless coupling constant a introduced above would 
be equal to a = (w/2)R/Rq, which can have values both 
larger and smaller than 1. 

Finally, for optical phonons, we use the Frohlich in- 
teraction Hamiltonian (Ref. 1831 1 with a dimensionless 
Frohlich coupling constant of a = 0.07 (GaAs) to obtain 

the rough estimate (FF) — 5{uj — ojlo) • (lmeV) 2 ■ 

(lOOnm/rfo), with ojlo ~ 5 • l0 13 Hz. This yields a z 
deviating from unity by about 10~ 3 . 

However, in the plots we have chosen z = 1/e for illus- 
trative purposes. 

The resulting behaviour of T((p, V) at T = 0, calcu- 
lated from Eq. l|T3jl . is shown in Figs. 0and[Sl In the 
case of constructive interference (ip = 0, Fig. QJ, the 



11 



decay rate for the "weak baths" (a,b,c) starts out from 
r/2r = z at V = and goes to T/2r = (1 + z 2 )/2 at 
eV/uj c 3> 1 . The initial deviation from the constant value 
of z at low voltages is given by the power-law l/ 2 ^ 1 ) con- 
tained in Eq. 1281 . In contrast, the decay rate for the 
Ohmic bath (d,e) starts at Y = 0, rising with a power- 
law and saturating at a value of r/2r = 1/2, corre- 
sponding to an equal admixture of odd and even states 
in the reduced density matrix of the electron coupled to 
the bath. For destructive interference (ip = n, Fig. EJ, 
the behaviour of (a) and (b) at low voltages is given by 
(see Eq. j28|l ). while the decay rate of the Ohmic 
bath (d,e) remains the same as that for tp = 0. In the 
special case (c) of the gapped bath, we observe perfect 
destructive interference up to the excitation threshold of 
the bath at eV = u> c , where r(<^ = w, V) increases in a 
stepwise manner for the first time, with the next increase 
at eV = 3uj c . Note that, on the other hand, Y(ip = 0, V) 
increases at even multiples of the excitation gap. The 
difference comes about because it is only the emission of 
an odd number of phonons into the bath that reveals the 
location of the electron, as discussed above. This feature 
would be absent if the two dots were coupled to two inde- 
pendent baths, whereas the other qualitative properties 
would remain the same. 



From the decay rates at <p = and ip — n, we may cal- 
culate the visibility v of the "interference pattern" that 
is defined by the dependence of Y on ip. The result is 
shown in Fig. 03 As we have noted before, the visibility 
is always zero for the Ohmic bath. On the other hand, 
for the "weak baths", it is perfect (equal to 1) at V — » 0, 
due to the perfect destructive interference, regardless of 
the suppression factor z appearing in Y(ip — 0). In gen- 
eral, the visibility decreases towards higher bias voltages 
before saturating at the limiting value of z 2 . However, 
in contrast to intuitive expectation, the decrease may be 
nonmonotonous, i.e. the visibility of the interference ef- 
fect may actually be enhanced by increasing the supply 
of energy available to the electron, although the decay 
rate Y always increases monotonously at any V. This 
is particularly striking in case (c), where the visibility 
drops down to zero in a certain range before rising again. 
The decrease down to the exact value of is related to 
the special choice of (^4> 2 ^ = 1 ( z = V e )i which gives 
equal strengths of the peak at E = and the first peak 
around E = u> c , which then are able to cancel in the 
integral 7- over P-{E) that is proportional to the visi- 
bility. However, the physical reason for a dip in visibility 
is rather generic: In that energy range, the decay rate Y 
for ip — 7r has already increased due to dephasing, while 
the blockade-type suppression of the value of Y for ip = 
has not yet been lifted. This is a consequence of the 
even-odd effect discussed above. 



VII. SEQUENTIAL TUNNELING THROUGH 
THE DOUBLE-DOT 

Up to now, we have discussed in detail the influence of 
the bath on the tunneling decay rate of an electron which 
has been placed onto the two dots in the symmetric su- 
perposition. In order to complete the picture, we have to 
calculate the sequential tunneling current through such a 
double-dot interference setup. This will be done by deriv- 
ing and solving a master equation for the reduced density 
matrix of the double-dot system, taking into account the 
system-bath coupling exactly, while the tunnel-coupling 
is treated in leading order. We are interested specifically 
in the nonlinear response, i.e. in how an increasing bias 
voltage helps to destroy the phase coherence. The tun- 
neling rates calculated previously will serve as input to 
the master equation. 

However, in order to facilitate the understanding of the 
results, we first turn to a qualitative description of the 
situation without the bath. 

At ip = 7r, tunneling is completely blocked, since the 
left reservoir only couples to the even state |e), while the 
right reservoir couples to the antisymmetric (odd) su- 
perposition, |o). At ip = 0, both reservoirs couple to |e), 
whereas |o) is completely decoupled from the leads (com- 
pare the discussion in Ref.|2g). This means that a current 
may flow if \o) is empty. However, if \o) is filled, the cur- 
rent vanishes, because double-occupancy is forbidden in 
our model. Since there is no way to change the occupa- 
tion of |o), the stationary density-matrix of the double- 
dot at ip = will be any convex combination of these two 
possibilities (at T — 0, in the absence of other relaxation 
paths). At any value of p in between these extremes, 
there is always the state |*) = (|+> - er lv \ -))/V2, 
whose decay into the right lead is blocked by destructive 
interference. As there is a nonvanishing overlap between 
|\&) and the state |e) which is reached by tunneling from 
the left lead, one will observe an accumulation of pop- 
ulation in until the current is blocked again. This 
argument holds at T = 0, while at finite temperatures 
the electron can decay towards the left lead and make a 
new attempt. Therefore, in this simple picture, the sta- 
tionary current at T = would be zero at any ip except 
for ip = 0, where it is undefined. 

However, one has to take into account that the cou- 
pling to the reservoirs does not only lead to decay but 
also to an effective tunnel coupling between |+) and |— ). 
Although this cannot change the blockade of the current 
at p = 7r (leading only to an energy shift of |e) vs. |o)), 
it does lift the blockade at other values of ip. This is be- 
cause the blocked state |^) is no longer stationary, such 
that an electron will not remain there forever. The de- 
generacy at ip- = still remains. Therefore, in the ideal 
case without coupling to a bath, we expect the current 
to vanish at ip = it and to rise towards a maximal ampli- 
tude near ip = 0. According to the previous argument, 
at T — this maximal amplitude will be determined by 
the effective tunnel-coupling between the dot states. 
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Introducing the bath will then lead to renormalization 
effects and spoil the perfect destructive interference at 
higher values of the bias voltage (or temperature) , quali- 
tatively in the same way as it has been explained above. 
We will show that the actual visibility vj of the current 
interference pattern I(ip) is given by a monotonous func- 
tion of the visibility v introduced above for the tunneling 
rate (at symmetric bias). 

We start with the Hamiltonian that is obtained after 
applying the unitary transformation of the independent 
boson model l|A5|l onto our Hamiltonian QJ: 

H' = e' (h+ +h-) + U'n+h-+H B +H L + H R + V' (36) 

Here e' is the (renormalized) energy of the two states, 
which we will take to be e' = from now on. U' is 
the interaction constant that involves both the Coulomb 
repulsion as well as the effective attractive interaction 
induced by the bath. We assume U, U' 3> T, eV, such 
that double-occupancy is forbidden. 

The term which we will treat as a perturbation is given 
by V', describing the tunneling to the left and the right 
leads in the presence of the bath. It is the transformed 
version of V (compare Eqs. <EJ| and {3} and AppendixEJ, 
where the additional fluctuating phase factors exp(±i<^>) 
have been introduced: 

v' = J2 Yl + h - c - > ( 3? ) 

j—l,r a— +,— 

where 

l ± = e ±l H (38) 
1 = (39) 

k 

f+ = e+% (40) 

f_ = e-^e itp f (41) 

r = (42) 

k 

As usual, the current through the device does not only 
depend on the rates for electrons to tunnel into and out of 
the dots, but also on the stationary state which the sys- 
tem assumes in the nonequilibrium situation, i.e. under 
an applied bias voltage. 

We will now derive a master equation for the reduced 
density matrix p of the double-dot system, which con- 
tains the populations p++,p--,poo ("0" denoting "no 

electron") and the coherences p^ and p |_ (with poo = 

l-p++-p — , p a o = Poa = Ofora ^ 0, andp_+ = p+_). 
We cannot simply use the standard kind of master equa- 
tion, since we have to deal with two degenerate levels 1+) 
and |— ), and it is important that a tunneling event may 



create a coherent superposition of |+) and |— ) (for exam- 
ple the even state |e)). Such a master equation - for de- 
generate levels - has also been employed in Ref . |n] (with- 
out coupling to the bath, and evaluated in the linear- 
response regime) . The equation is different from that em- 
ployed in the "orthodox" theory of sequential tunneling, 
where no coherent superpositions are involved. Note that 
for a finite tunnel-coupling the levels could be treated 
as degenerate as long as their energetic distance is much 
smaller than the level-broadening due to tunneling. How- 
ever, as we consider the limit T — > 0, we need to have 
exactly equal energies. Otherwise, the energy of the hole 
that is created in the left electrode would betray the dot 
state which the electron has entered, thus preventing any 
coherent superposition to form. 

Given the initial reduced density matrix p(0), and as- 
suming the state of the environment (bath and reservoirs) 
to be independent of the electronic state on the dot at 
t = 0, we obtain the time-evolution p{t) by tracing over 
the environmental degrees of freedom ("E"): 

= m 

/ dh / dt2tr E [V'(ti)V'(t 2 )p(0)®pE + h.c.} 
Jo Jo 

+ f dh f dt 2 tr E [V'(h)p(0)®pEV'(h)] 4(43) 
Jo Jo 

Physically, by using factorized initial conditions, we 
neglect correlations between subsequent tunneling events 
which could be due to excitations in the electrodes or in 
the bath: Since the tunneling rate is very small, these 
excitations will have traveled away from the double-dot 
until the next event takes place. The entanglement be- 
tween electron and bath (discussed in the previous sec- 
tions) would preclude factorized initial conditions, if it 
were not treated indirectly in this approach (via the uni- 
tary transformation) . Note that we do not have to make 
any secular approximation at this point, unlike the usual 
derivation of a master equation^!. It turns out that all 
contributions only depend on the time-difference t\ — t 2 
anyway, because the dot levels are degenerate. There- 
fore, in the long-time limit t — > oo, the integration over 
(ii + t2)/2 results in a factor t, and the endpoints of the 
integrals over t\ — ti may be extended to oo. This yields 
the desired master equation that will determine the sta- 
tionary p, as well as the current, in the limit of weak 
tunnel coupling. 

In the expectation values of products V'V' only those 
contributions remain which combine d a j a (tunneling out 
of the dots) with jtdt (tunneling onto the dots): 




+dl(s)d P p(]i(s)j p ) + h.c.} 
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+ 00 



ds [d a {s)pdl (jp3a(s)) + 

dUs)pd (j jU S ))} . (44) 



(Note that there is no minus sign from fermion opera- 
tor re-ordering in this factorization of dot and reservoir 
part, as the reservoir fermion operators are dragged past 
an even number of dot operators; compare e.g.'®; alter- 
natively, it is also possible to define them as commuting 
operators, since there is no interaction between them). 
We get for the individual matrix elements (for brevity, 
the summation over j = l,r is implied): 



-P++ I ds {j f + (s)j+ 

ds (j+j+(s) 
it 



— oo 
r+oc 

+POQ 

' — oo 

fOO 

f I ds (j J + {s)j-) - h.c. 



/'+ = -/'- - J ds {j[(s)j+ 



(45) 



ds (f-j-(s] 



-/'no / ds Cj-3+(s) 

I — oo 

-P++ J ds (j\j-{s) 



-P— 



ds (j+(s)j- 



(46) 



The equation for p__ follows from that for p++ by 
interchanging indices + and — . 

Now we have to evaluate environment correlators, such 
as the prefactor of p++ in the second equation (e.g. for 

j = r): 



(47) 



By introducing the bare tunneling rates T R / L ^ 



2ttD 



ML) 



(compare Eq. ©), we get, using 



Eq. 1421) (remember r creates a reservoir electron): 



^f(s)) = ^° j de (1 - Me)) e+*" . (48) 

Here we have neglected any energy-dependence of the 
tunnel-coupling and electrode DOS, assuming the rele- 
vant voltages and temperatures to be sufficiently small 
(but see below) . The bath correlator in I|47J1 evaluates to 



cxp(— yj)(j)(s)j — (cjrp, which can be expressed by us- 
ing the definition lfr2)l for P_(w). There, we have to set 
-s because of the reversed order in the ^-correlator: 



(49) 



Therefore, we obtain: 



ds /f+f- (s)y 



e ^lM j de( i_ /fl(e))F * ( _ e)j 



(50) 



with 



P_(e) = iy dujP-{Lu) J dse l{t - u)s = 

P_(e) + - / 

7T J e-uj 



(51) 



The integral in the second line is understood as a 
principal-value integral. In order to abbreviate expres- 
sions like this, we introduce the following definitions for 
the effective in- and out-tunneling rates as well as the 
effective tunnel couplings generated by the electrodes: 



IL(-) 


-!■„/ 


iL(-) 


-r„/ 


A L 


= _£i» 




7T 


It- 


= 7L-[P 



p_H 

6 + U) 



iA L . 



(52) 
(53) 

(54) 
(55) 



Analogous definitions hold for L i— > P. Eq. l|52)l is 
equivalent to the definition itHjl used for 7(_) in previous 
sections. Note that the effective tunnel coupling 
depends on P_ , because it arises from transitions between 
the states |+) and |— ), via an intermediate lead state. In 
the expression for A £ ( H v the energy-dependence of the 
density of states and the tunnel coupling to the reservoir 
electrode should be kept in order to have a convergent 
integral. We will take this into account by introducing 
an effective upper energy cutoff A in the integral. Using 
these definitions, 15011 is equal to exp(i(p)j^ i _/2. 

One might wonder why the effective tunnel couplings 
A_L(i?) do depend on the occupation of electron states in 
the reservoirs. After all, in the non-interacting case, it is 
possible to calculate such a change of the effective single- 
particle Hamiltonian prior to filling in the electron states. 
Alternatively, in a calculation that already takes into ac- 
count occupation factors, there would be two contribu- 
tions which add up to an integral that does not depend 
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on the Fermi function. However, we consider the inter- 
acting case U = oo, such that (even without the bath) 
one of these contributions is missing (since it would in- 
volve intermediate states with double occupancy). The 
resulting logarithm is analogous to that which appears 
in the Kondo problem. This effective tunnel coupling 
has also been discussed in Ref. H3, for the case without 
a bath. There, the upper cutoff A was provided by the 
Coulomb coupling U, since for higher energies double- 
occupancy is no longer forbidden and the non-interacting 
case takes over (where two contributions arise that can- 
cel each other). If we take the limit U — > oo, then A will 
be set by a cutoff in the tunnel matrix elements (or the 
electron reservoir's density of states). 

The general master equation for the reduced density 
matrix of the double-dot, derived in the limit of weak 
tunnel coupling but arbitrary electron-bath coupling, fol- 
lows by inserting the definitions l|52l) - l|55)) into Eqs. jUt 
and 

P++ = -p++(lL+Jn.) 

+Poo(lT+lh n ) 

-^{e^ R -+%-)-h.c, (56) 

P— = -P—(il+Jr) 
+Poo(7f +7k") 

-^(e-^R-+%-)-h.c, (57) 



where p^ = /?++ = p Taking into account that we 

are working in a transformed basis, this describes just the 
entangled state , confirming the starting point of our 
earlier discussion. Note that the out-tunneling contribu- 
tion for p++ also depends on p_| , for example. This 

reflects the fact that a superposition between the two 
states may be blocked from decaying into the lead, while 
each state separately can decay. 

The stationary density matrix is obtained by demand- 
ing dp/dt = (and using the relations poo = 1 — P++ — 
P-- and p |_ = p* ). This will give us the density ma- 
trix in zeroth order Tq in the bare tunnel coupling, which 
we need to calculate the current in leading order Tj. 

We can obtain the current from the contribution of the 

left electrode to the change p ++ + p in the double-dot 

occupation (i.e. keeping only terms that stem from the 
left electrode in the master equation). This is equal to 
the right-going current in the stationary limit: 



- = (f>++ +P—)l = 
e 

2poo7i" " 1l(p++ + P-) - 2 7i _i?e[p+_] . (60) 

An alternative way of deriving the current would be to 
start from the general Meir-Wingreen formula^ which 
expresses the current in terms of the exact Green's func- 
tions of the double dot, to be calculated in presence of 
the tunnel-cou plin g and the bath. This has been the ap- 
proach of Ref. O for the case without the bath, and we 
have checked IfiOt to give the same result in that case. 



P+- = -P+-{1L +1r) 

i / iip in , in \ 

-^(e*7*-+7i-). (58) 

The ingredients of the master equation obtained here 
may be interpreted as follows: 

One part of the right hand side corresponds to the uni- 
tary time-evolution generated by the effective tunneling 
Hamiltonian 

Hj ff = i(e*A„ + A L ) |+) (-| + h.c. . (59) 

Furthermore, the in-tunneling contributions in the 
equations for p++ and depend on P(E), while that 

for p^ is determined by P-(E), since it describes the 

creation of a coherent superposition of |+) and |— ) (which 
is hindered by the bath). This term would be absent in 
the usual master equation. In particular, if j™_ ~ * 7™> 
which will be the case at T = for vanishing bias between 
the dots and the left electrode, an electron tunneling from 
the left lead will end up in the coherent superposition 



VIII. EVALUATION OF THE SEQUENTIAL 
TUNNELING CURRENT AND THE VISIBILITY 

In order to evaluate the current as a function of tem- 
perature T, bias voltage V and phase difference <p, we 
will now specialize to the case of symmetric bias and left- 
right symmetric tunnel couplings (r^o = ^lq = To). All 
essential features (in particular the perfect destructive in- 
terference in absence of the bath) are independent of this 
assumption. We will find that the current is symmetric 
under ip i— > —ip even for the nonlinear response consid- 
ered here, due to the symmetry of the model (compare 
Ref. O for a systematic analysis of phase-locking in a 
variety of interference geometries) . 

We find from Eqs. JSJ-lO, usin & /( e ) = 1 ~ f(~ e ) : 

7fl(-) = llU = 7(-) = To J de f(e - p) P ( _) (e) , (61) 

where p — eV/2 is the chemical potential of the 
left reservoir. This is definition 11411 . with eV replaced 
by p — eV/2 (since we deal with the symmetric bias 
case). Furthermore, we use the condition of detailed bal- 
ance, P(_)(— E) = exp(—f3E)P(_* ) (E) (see, for example, 
Ref. 1811 1 , which leads to 
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Figure 7: The current I for different values of the visibility 
v = 7-/7 = 0.8, 0.9, 0.99, 0.999, 0.9999 (from top to bot- 
tom). The limits <p — > and v — > 1 do not commute. Other 
parameters held fixed: A = e - ^ = 0.2 and 5l = 8r = — 1. 



7l(-) = 7ii(-) 



e-^ 7 (-) 



(62) 



The effective tunnel couplings are still different (be- 
cause of the different Fermi distributions) : 



■--/ def(-(e T n)) [du 



e + w 



(63) 



The lower sign belongs to the right electrode. 

For the special case of T — 0, electrons always enter 
from the left and go to the right, such that we have = 
1l- = 7r = 7h_ = and Tr(-) = 7i"_) = 7(-) s with 



(64) 



7(-)=r / deP ( _)(e). 
Jo 

The effective tunnel couplings are, at T = 0: 



A 



L(R) 



— J du)P-(u)) In 



A + w 
\fi±u\ 



(65) 



Note that, without any bath present, &l(r) will have 
a logarithmic singularity at /i — > 0, for T = 0. The upper 
cutoff A will be given by the minimum of the Coulomb 
repulsion energy U and the bandwidth of the reservoir's 
electronic energy band (or by some cutoff in the tunnel 
matrix elements). For the purposes of our discussion, we 
assume A ^> /i, ui. 

In the the limit of high bias voltages (u <C A, fx), we 
obtain effective tunnel couplings whose magnitude goes 
as z 2 and decreases logarithmically with increasing li: 



Ar « A? 



---In 



-1 / dcjP_( W ) = -z 2 ^ln[- 

M i t La* 

(66) 

By solving the master equation for the stationary den- 
sity matrix and inserting the result into Eq. Il(i0l) . we 
obtain the expression for the current through the double 
dot in terms of all of the quantities mentioned previ- 
ously. In general (at arbitrary T), it is found that the 
current may be written as the product of 7 with a di- 
mensionless function of the phase difference ip and the 
ratios v = 7 _/ 7 , S L{R) = A z ,( fl) / 7 and 0^i: 



(67) 



The complete expression for I is very cumbersome, 
although it may be found analytically by straightforward 
solution of the master equation (it is listed for T = in 
AppendixEJ. Therefore, let us first discuss the situation 
without coupling to a bath. In that case, we obtain 



S L = S R ^S = -^ ( A -/(A*-e) (68) 

7f J-oo e 

and 7 = 7-= Tof{— fi). The current turns out to be 
(with A = er^): 



I 



4(1 



A)Q5 2 + A)cos 2 (f) 



3<5 2 + 2(1 + A + A 2 ) + 36 2 cos(v?) 



(69) 



Several points should be noticed about this expression: 
Firstly, the destructive interference at 99 = it remains 
perfect regardless of temperature, because there are no 
current-carrying states at all. At zero temperature (A = 
0), the maximal amplitude of the current is I m ax/ej = 
25 2 /(35 2 + 1), which vanishes when the effective tunnel 
coupling S goes to zero. This has been explained above 
as a consequence of the possible transition into a current- 
blocking state, which can only be undone by the effective 
tunnel coupling. At finite temperatures (A > 0), the 
maximal current is nonzero even for S — > 0, where it 
approaches the value oi I ma x/ej = 2A(1 — A)/(l + A + A 2 ). 
This has a maximum at around T ~ fx. It vanishes for 
larger temperatures as fi/T, which is to be expected for 
tunneling through a localized level (decreasing derivative 
of the Fermi function). In addition, the shape of I(<p) 
depends on S and A, with a sharper minimum at <p = ir 
in the case of larger \S\. In the limit of 5 — > 0, the current 
becomes a pure cosine. At finite temperatures (as well as 
for v ^ 1) the behaviour is similar, except for the finite 
amplitude of the current at 8 — > 0. 

Now we turn to the situation including the bath. The 
general expression for the current is very lengthy, and we 
will omit it here. However, it turns out that the maximal 
and minimal current are functions merely of v and A = 
e - ^, while they are independent of 5l,r- 
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The amplitude of the minimal current (at tp 
given by 



I(tp = tt) 2(l + A)(l-A 2 )(l-i> 2 ) 
e 7 _ 3(1 + A) 2 + (1-A) 2 w 2 

while the maximal current (at tp = 0) is 



7T IS 



(70) 



l(p = 0) 



1(1 -A) 



(71) 



It should be noted that the expression l(69jl for the current 
in the ideal case seems to contradict this simple formula. 
However, that is because the limits tp — > and v — > 1 do 
not commute. This is shown in Fig. It means that 
for T = and 5l,r — > the maximal current calculated 
according to j71|. which is independent of <5l,.r, and the 
"typical" amplitude of the current (oc <5 2 ) may deviate 
strongly. The peculiar behaviour near tp = seems to 
be connected to the physical degeneracy of the case tp = 
0, v = 1 which has been discussed above. 

From these formulas, we obtain the visibility, defined 
in terms of the current: 



0.8 



0.6 



0.4 



0.2 



-- optical phonons 
— acoustic phonons 




eV/2w 



Figure 8: The visibility vi of the pattern I(<p), for piezoelec- 
tric coupling to acoustic phonons (b) (solid line) and for the 
optical phonon bath (c) (dashed line), plotted vs. fi = eV/2, 
at different temperatures T/lu c = 0.01, 0.05, 0.1, 0.2, 0.4, 0.5 
(top to bottom). Inset depicts energy diagram for tunneling 
in this situation. 



Vi 



I(tp = 0) - I(p = tt) 
I(tp = 0)+I{p = tt) 



(72) 



It can be expressed entirely by the visibility v defined 
previously in terms of the tunneling rates (Eqs. iflfil) . 

as well as the temperature-dependent factor A = 
e"^ (fi = eV/2): 



vi 



2(1 + A + A> 2 



3(1 + A) 2 - (1 + 4A + A 2 )u 2 



(73) 



This is a monotonous mapping of v to the interval 
[0,1], with only a weak dependence on A. The other 
parameters 8l,5r only modify the amplitude and shape 
of the current pattern I(tp). Therefore, all the statements 
about the visibility made in the previous discussion of 
the tunneling decay out of the symmetric superposition 
continue to hold up to this monotonous transformation 
(and with eV replaced by /i = eV/2). In particular, at 
T = 0, we have 



vi 



2v 2 



3 - v 2 



(74) 



The dependence of the visibility vi on the bias voltage 
eV = 2/i, the temperature T and the bath spectrum is 
displayed in Fig. |H1 for bath spectra of type (b) and (c). 
The decrease of Vi at /x = with increasing tempera- 
ture T in case (b) is well approximated by Eq. l(35|l for 
v(T, V — > 0) (employing the relation vj — v 2 /(2 — v 2 ) for 
fi = 0). (The functions Pf_-)(E) for finite temperatures 



have been calculated numerically using the fast Fourier 
transform, from the defining equation 11121) ). 

Note that for bath spectra with z — (i.e. exponent 
s < 1 at T = and s < 2 at T > 0) the visibility van- 
ishes entirely (at any V), as has been explained in the 
previous sections. We have already pointed out that this 
picture is expected to change if one treats the tunnel- 
coupling to higher order. However, we have to leave this 
analysis for the future. One possible approach to a non- 
perturbative (but still approximate) treatment of both 
the tunnel-coupling and the system-bath coupling at the 
same time seems to be the numerical "real-time renor- 
malization group" schema^. 



IX. CONCLUSIONS 

We have analyzed dephasing in tunneling through two 
parallel single-level quantum dots with a fluctuating en- 
ergy difference between the dots. The disappearance of 
perfect destructive interference in a symmetric setup has 
been taken as a criterion for "genuine" dephasing, as op- 
posed to mere renormalization. The coupling to the bath 
has been taken into account exactly, via the "independent 
boson model" and the concepts of the U P(E) theory" of 
tunneling in a dissipative environment, while the tunnel 
coupling has been treated in leading order. 

We have discussed in detail the behaviour of the den- 
sity matrix of a single electron that has been placed in 
a superposition of the two dot levels. The bath mea- 
sures (to some extent) the position of the electron, such 
that the electron's density matrix becomes mixed. How- 
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ever, this allows direct conclusions about the "incoherent 
current" only in the limit of high bias voltages, corre- 
sponding to a fast "projection" measurement of the elec- 
tron's state. For lower voltages, only the low-frequency 
part of the bath spectrum contributes to the lifting of de- 
structive interference. Thus, for any "weak bath", whose 
spectrum falls off fast towards low frequencies, the vis- 
ibility of the interference effect becomes perfect in the 
limit of low bias voltages V and temperatures T, when 
the energy supplied to the electron is vanishingly small. 
This is the case for a fluctuation spectrum oc oj s with 
s > 1 (s > 2) for T = (T > 0). The visibility may 
show a nonmonotonous behaviour as a function of bias 
voltage. For "stronger" spectra (smaller exponent s), in- 
cluding the Ohmic bath (s = 1), there is the well-known 
zero-bias anomaly (suppression of the tunneling current 
at low voltages), which affects equally both the cases of 
constructive and destructive interference. Therefore, the 
visibility vanishes exactly at any bias voltage in our ap- 
proach, where the tunnel coupling has been treated only 
in leading order. Although there is always a suppres- 
sion of the magnitude of the tunnel current for the case 
of constructive interference, this may be interpreted as 
a mere renormalization of the effective tunnel-coupling, 
since the perfect destructive interference is not affected 
and since it occurs even for a bath with an excitation gap. 
The full dependence of the sequential tunneling current 
I(ip) on voltage, temperature, bath spectrum and phase 
difference <p between the interfering paths has been de- 
rived by setting up a master equation for the state of the 
double-dot (which is special due to the degeneracy of dot 
levels) . 

The major questions that have remained open in our 
analysis are related to the behaviour at stronger tunnel 
coupling. In particular, the perfect destructive interfer- 
ence may also be overcome by correlated tunneling of sev- 
eral particles (with an intermediate "virtual" excitation 
of the bath), and this process will therefore contribute 
to dephasing, although it is expected to be suppressed 
strongly at low voltages and temperatures. Likewise, 
the visibility for the Ohmic bath (or other strong baths), 
which turns out to be zero in the present approximation, 
may be changed at low bias voltages and temperatures 
comparable to the tunneling rate. This will require other 
methods to analyze the competition between strong tun- 
nel coupling and system-bath coupling. 



Appendix A: INDEPENDENT BOSON MODEL 

For reference purposes, we describe here the canoni- 
cal transformation employed in the independent boson 
model. See RefUS for more details (concerning the case 
of at most a single particle). Consider a set of electronic 
levels j that couple to bath operators Fj which are as- 
sumed to be linear in the coordinates (and momenta) of 
a bath of harmonic oscillators, Hb- 



(Al) 



did., 



Here £j is the unperturbed level energy and hj — u, .^ 
is the number of particles on level j. The fluctuating 
fields are characterized completely by their power spectra 
at T = 0, 



T=a i 
u 2ir 



+ °° I „ „ \ T=0 

dte 1 ^ {Ft^Fj 



(A2) 



Here we will restrict ourselves to the case where the dif- 
ferent variables commute, [F;, Fj] =0. As a consequence, 

the spectrum (FiFj) is real- valued, but there may 

still be correlations. 

The most straightforward solution proceeds via a uni- 
tary transformation^ (essentially a gauge transforma- 
tion). One introduces the fluctuating phases (f>j, whose 
time-derivatives are given by the Fj : 



i[H B ,<Pj] = -Fj 



(A3) 



The exponent generating the unitary transformation is 
defined as: 



(A4) 



Applying the transformation to the Hamiltonian in Eq. 
llA"T)l yields: 



H' = e~ lx 



He +i * = £ j fl j - Jlfi&i + H b ■ (AS) 



ij 
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The coupling between system and bath has been elim- 
inated, resulting in an effective interaction between par- 
ticles on the different levels, with: 



Ju 



duj 



T=0 



(A6) 



The Jij are real-valued and independent of tem- 
perature. For I = j they describe energy shifts of 
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single-particle levels. The canonical transformation also 
changes the particle annihilation and creation operators, 



bath. 



it, 



3 ' 



(A7) 



This will affect all Green's func- 



and d'J „ r 

tions and, therefore, also the time-evolution of the single- 
particle density matrix. In addition, it becomes impor- 
tant if a tunneling part is added to the Hamiltonian, 
where the operators d*- appear, such that they have to 
be transformed according to l|A7l) . However, since the 
phases 4>j and the particle operators commute (even 

at different times, when evolved according to H'), the 
evaluation of Green's functions always splits into a part 
referring to the particles and a separate average over the 
bath operators. This is the major simplification brought 
about by the "diagonal coupling" between system and 



Appendix B: CURRENT EXPRESSION FOR 
SEQUENTIAL TUNNELING THROUGH THE 
DOUBLE-DOT 



At T = 0, for the symmetric situation, the current / is 
riven by / = ej I [v, 5 L , Sr], with: 



Io[v,S l ,Sr] = 2-[ 



-6 2 +(v 2 -l)(l + S 2 R ) + 
28l$r(v 2 — 1) cos ip + 5\v 2 cos 2 ip] ■ 
[-3(5| + 2S L 5 R v 2 + (1 + 8 2 R )(v 2 - 3) + 
2{v 2 {\ + S 2 L + 8 2 R ) + 6 L S R (v 2 - 3)) coscp + 



S L (S L + 2S R )v 2 cos 2 v?]" 



(Bl) 
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